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Abstract 



We solve the Schwinger-Dyson equations for (2+l)-dimensional QED in the presence of 
a strong external magnetic field. The calculation is done at finite temperature and the fermionic 
self energy is not supposed to be momentum-independent, which is the usual simplification in 
such calculations. The phase diagram in the Temperature-Magnetic field plane is determined. 
For each value of the magnetic field the critical temperature is higher than in the constant mass 
approximation. In addition, the latter approximation shows a relative B-independence of the 
critical temperature, which occurs for stronger magnetic fields in the momentum dependent 
case. 
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1 Introduction 



The mechanism of dynamical mass generation in the presence of a magnetic field is particularly 
important in connection with the vacuum structure of non-abelian field theories such as QCD; 
also the QED case is interesting both as a theoretical "laboratory" and because of possible ap- 
plications. In particular, scenaria of dynamical gauge symmetry breaking in three-dimensional 
QED U lead to interesting and unconventional superconducting properties of the theory after 
coupling to electromagnetism ||, and therefore may be of interest to the condensed matter 
physics, especially in connection with the high-T c superconductors. 

The phenomenon of magnetic catalysis has been studied by several groups || [|, |] in 
various models. It has been found that a homogeneous magnetic field induces dynamical mass 
generation even for the weakest attractive interaction between the fermions. In addition, in 
2+1 dimensions the mass generation is not restricted to a small number of flavours (which is 
the case in the absence of the magnetic field). 

The treatment of the problem necessarily involves simplifications; a serious one is the 
constant mass approximation, according to which the fermionic self-energy is supposed to be 
momentum-independent. There has already been a first attempt [EJ to treat the problem 
of dynamical mass generation for QED in 2+1 and 3+1 dimensions at T = 0, taking into 
account the momentum dependence of the fermion self-energy. The results have shown that 
there are differences (very important ones in 3+1 dimensions) showing up when the momentum 
dependence is taken into account, thus establishing the necessity to go beyond the constant 
mass approximation. In that work, attention was restricted to the strong magnetic field regime, 
to facilitate the calculations. In the present work we will extend the calculations in || to the 
finite temperature case in 2+1 dimensions. Extensions to thermal 3+1 dimensional QED will be 
deferred to a forthcoming publication. In the present work, only the regime of strong magnetic 
fields will be studied, to render the problem tractable. 



2 Fermions in a constant magnetic field 

To fix our notations we shortly review here the characteristics of fermions in a constant external 
magnetic field in 2+1 dimensions at zero temperature. The model we are going to consider is 
described by the Lagrangian density: 

£ = --F^F^ + i^D^^/ - m*tf, (1) 

where = <9 M + iga^ + ieA^ 1 , a M is an abelian quantum gauge field, F^ u is the corresponding 
field strength, and Aff 1 describes possible external fields; in this work A e * 1 will represent a 
constant homogeneous external magnetic field. Notice that the fermions feel both the quantum 
and the external gauge fields, however we have allowed for different coupling constants, g and 
e, in order to give an effective description of condensed matter systems 0. We recall the usual 
definition g 2 = Ana. 

We will choose the "symmetric" gauge for the external field 

A^(x)=0, = -~x 2) Af\x) = +^x 1 (2) 
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for which we know from the work of Schwinger |7j that the fermion propagator is given in the 
Minkowski space by the expression: 



S(x,y) 



^iex^ Af i xt (y) 



S(x - y), 



(3) 



where the translational invariant propagator S has the following Fourier transform: 



Sip) 



2 „2 tan(|eB|s) 



X 



(p°7° + m)(l + TT tan(|e£|s)) - p^il + tan\\eB\s)) 



(4) 



where p L = (p 1 , p 2 ) is the transverse momentum and a similar notation holds for the 7 matrices. 

Let us now turn to the finite temperature case. We will denote the fermionic Matsubara 
frequencies by vbi = (21 + 1)ttT and the bosonic ones by wi = 2lnT. The translational invariant 
part of the bare fermion propagator can be expressed in Euclidean space by performing the 
rotations po — > iwi and s — > —is/\eB\: 



x 



— + m)(l — «7 1 7 2 tanhs) — p J "7" L (l — tanh 2 s) (5) 



Si(p±) = - r - B r / dse 
\eB\ Jo 



The Euclidean 7 matrices satisfy the anticommutation relations {7^, 7^} = —25^, with fi, v = 
1,2,3. 

The fermion propagator has another representation in the lowest Landau level (LLL) 
approximation ((§f,0), which at finite temperature reads: 

§i(p±) = -ie-^/l eB l "fj 7 ' + 2 m (l - i^signieB)) . (6) 

Inspired by (|6]) and taking into account the fact that the dressed fermion propagator has the 
same phase dependence on A ext as the bare one 0: 

G(x,y) = e iex " A » t{y) G(x-y), (7) 

we make the following ansatz for the Fourier transform of the translational invariant part of 
the full propagator: 

G,(px) = -* A ' m Z §ff^f (1 " nVsign(eB)) (8) 

where Z\ is the wave function renormalization and Mi the dynamical mass of the fermion. 
Both quantities depend only on the Matsubara index; this is a consequence of the restriction 
to the strong field regime, where the LLL approximation reduces the fermion dynamics to 1- 
dimensional dynamics ||. We will study the dynamical generation of fermionic mass only in 
the case m — 0. 
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3 Integral equations and the recursion formula 



We can write the equations satisfied by the wave function renormalization Z n and the dimen- 

M n j \J\eB\ from the corresponding relations obtained from the 
(|PH|, 0) making the substitutions 



sionless dynamical mass /i r 
Schwinger-Dyson equation at T 



Vz 



wi and 



dp3 
2vr 



TJ2 

i 



to obtain 



2at ^ ZiCji 

Z„ = 1- 



n i Zfuf + fif Jo 

= 2dt^ _ 2 ^' 2 f°° rdre- r2 ' 2 V n ^{r) (9) 

where we introduced the notations a = a./ J\eB\, t = T/\J\eB\, &i = Wi/y\eB~\ and r stands 
for the dimensionless modulus of the transverse momentum: r 2 = qj_/\eB\. T> m (r) is the 
dimensionless longitudinal component (//, z/ = 3, 3) of the photon propagator D^ u which is the 
only one to play a role since in the LLL approximation the fermion propagator is proportional 
to the projector operator P = (1 — i r j lr j 2 sg{eB))/2 and the original Schwinger-Dyson equation 
contains the product 

P-fP-fPD^ = -PD 33 (10) 

We note that for I = n in @ the integration over r is divergent in the infrared if we use the 
bare photon propagator. Therefore we will use the dressed photon propagator, which is the 
subject of the next section. In figure we show the quantity (r 2 + u 2 )V n (r) versus r. For the 
bare photon propagator this quantity is 1; in the figure we have set n = (for which the bare 
photon propagator leads to a divergence). We observe that the dressed propagator equals the 
bare one for most of the range of r (from about 3 to infinity). The difference of the two is 
restricted to a small neighborhood of r = 0. This is how this propagator cures the infrared 
divergencies associated with the use of the bare photon propagator. 

We will calculate the longitudinal component of the polarization tensor using the full 
Schwinger representation (|J) of the fermion propagator since in the LLL approximation it 
is zero, both at T = and T > 0. This computation will be valid for any magnitude of 
the magnetic field and will thus include the strong field approximation. We only need the 
longitudinal component of the polarization tensor since we have [IT| (with the dimensional 
reduction from 3+1 to 2+1) 

VJr) = - n - n =P r 33 + , „ ^ ^ (11) 



r 



2 + u 2 -F L (r 2 + cu 2 ) 2 



with 



?? = -^- 2 and F=^l where fg(r) = ggM^ , (12) 
r z + PI \eB\ 
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which leads to 

1 r 2 fl 33 (r) 
X> ( r ) = -| n 1 ; - (13) 

r 2 + ui 2 n ( r 2 + CJ 2) 2 ( r 2_n33( r )) 

in the Feynman gauge that will be used in this paper. 

To solve (^|), we will proceed as in the case T = 0: we put a trial series ^ %al in the 
integral equations and obtain a first estimate fj,^' of the series fx n . We then use this series in 
the integral equation and build up a (converging) iterative procedure. The trial series will be 
a solution of the recursion formula that we derive now, following a procedure similar to the 
one that gave the differential equation at zero temperature in ||, 0. We start by splitting the 
summation over / and write, starting from equation @ satisfied by /i n : 

^ ~ 2 ^E V 2^\ 2 f°° rdre- r2 / 2 V n (r) 

\l\<n Z l U l + J ° 

+ 2» £ 7V ,?' ,2 C rdre^l^r) (14) 
such that the difference fi n+ i — fi n reads 



Hn+l - fi n = 2atf n ^ v2 X\ 2 ( 15 ) 

\l\<n L l U l + 



where we defined 



roc „ 

/„=/ rdre- r ' 2 [V n+1 {r) - V n {r)} (16) 
J o 
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Figure 2: ji n versus n for a = .01 and t = .001 



Since we wish to get rid of the sum in equation fllSQ , we eliminate it manipulating the expressions 
for the differences (i n +i — \i n and /i n+ 2 — AWi- We finally obtain: 

7 If — T /in+1 + T = 7 2 r, 2 + ,/ 2 (17) 

Jn+1 \Jn+l In) Jn Z/ n+1 Co> ra+1 -+- fl n+1 



The solution of ( fPTD is found by giving (equal) initial values to /xo an d A*i, such that lim n ^oo fj, n 
0. We plot in figure || the solution of ( |17D (where we take Z n — 1 for every n) as well as the first 
iterations in the integral equations (|9|). We see that the convegence of the iterative procedure 
is very quick. We show in figure |3| the wave function renormalization versus n, for which the 
convergence is also very quick. From the technical point of view, the function f n has been 
computed with a Gauss-Hermite quadrature of order 40 and the series /i n and Z n have been 
truncated to an order depending on the temperature, adjusted in such a way that a given 
precision was kept throughout the numerical computations. 



4 Longitudinal polarization tensor 

We now proceed with the calculation of the polarization tensor using the Schwinger represen- 
tation of the fermion propagator fl5|). 

The one-loop polarization tensor is 

IC(*l) = -47mT£ J ^tr{ 7 ^(^)7^-n(p± - k ± )} (18) 

We note that the A^ 1 - dependent phase of the fermion propagator does not contribute to the 
polarization tensor since in coordinate space this phase contribution is 
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Figure 3: Z n versus n for at = .01 and t = .001 



exp {ie (x^Afiy) + y^Afix)) } = 1 (19) 

as can be seen from the potential (|J). We also remark that, since we use the full expression for 
the fermion propagator, the result of this calculation will be valid for any value of the external 
magnetic field. 

Using the expression (|5|) of the fermion propagator, we obtain for the longitudinal com- 
ponent, after the integration over the transverse momentum p± 



-4aT 
\eB\ 



E 



dsda 



o tanh s + tanh a 



tanh s tanh a 



x 



+ ( 



tanh s tanh a 



(tanh s + tanh a) 2 ytanhs + tanh a 
■ 1 ' ' 1 + tanh s tanh a) 



tanh S tanh g _ 1 r/ + ) /^2 + 2) + / Wn -2w ;)] 

tanh s + tanh cr |eB| LV ^ v ( ' « t/J 

e.B| \ 

^1 — tanhs)(l — tanha) 



m 



(20) 



where we will take the fermionic dynamical mass for m and consider that it does not depend on 
the momentum. We use this approximation since it is the only way to keep the calculation of the 
polarization tensor tractable. Without this assumption the momentum integrals cannot be done 
analytically; the numerical treatment on the other hand would be prohibitively demanding. 

As in [fL2|, we make the change of variable s = u(l — v)/2 and a = u(l + v)/2 to obtain 



-2aT 
\eB\ 



E 

i 



udu I dve 1« 



J_ cosh u — cosh uv 



-j^[tf+rn 2 +(l+v)w n {w n /2-Cj l )\ 



-1 
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m ) coth u — 



\eB\ 



sinh u 



cosh u — cosh uv 
2 sinh 3 u 



:2V. 



We integrate by parts: 



due~^ u) - 



u 



sinh u 



Jo \ du J 



(22) 



where we have discarded the surface term . This term is finite but would lead to an infinite 
summation over the Matsubara modes. We then obtain the final expression: 



-2a 



1 r.2 

' J_ cosh it — cosh uv u r 2 j_ l — v z 2 



du I dve • /; " ■ i: - 



2 

4 ^nj 



(23) 



x T £ e -^ 2 



/c 2 cosh uv — v coth u sinh ut> 



sinh u 



+ cothu (2«W 2 + uvw n Wi - |efi|) 



where W 7 / = — - 2 V w n . 

At this point we must make an important remark. In equation (|23"D there is a potential 
divergence of the integral over u coming from the sum involving \eB \ cothw. We find out that 
if we perform the summation over the Matsubara modes before performing the integration, this 
would-be divergence cancels against the sum involving 2iiW^ 2 cothii. This can be easily seen 
using the Poisson resummation [ |T3|j : 

OO / \ 1/2 OO 22 

J2 e~ a ^ 2 = (-J J2 e~^- 2mz \ (24) 

£ = — OO l = — OO 

which shows that the difference below has no divergence since 



T £ (2«W? - |eB|) e-iafr = ^ ' ^(-l^/V^ cos[vrn/(l + v)\ (25) 

i=-oo ^V 7 ™ 1 i>i 

so that the integration over u in ([23|) is safe, both on the infrared and the ultraviolet sides. 
The conclusion is that one should first sum over the Matsubara modes and perform the integral 
over u afterwards. 

If we take the limit T — > in (^3|) , we recover the zero-temperature results given in [|J 
since the substitutions Wi — > ps and TJ2i^ (2tt) -1 / dp 3 lead to 




lim TYuvw n Wie-^\ w ? = (26) 

such that (c<j n — > k 3 ) 
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„„ , , — 01 f°° , fc x cosh u cosh uv u [„i i 1—V? u3] 

\imUf(k ± ) = — . / duy/u dve ^ 2sinh " TW m + 4 fc 3 i 

T^O l/dpRMO J-i 



2Jn\eB\ J o 



, 9 cosh uv — v coth u sinh uv 

x k 2 ± — (27) 

sinhw 

We call the attention of the reader to a rather tricky aspect of these limiting procedures. In 
equations (|26|) we have taken the limit T — ► before we integrate over w and the result has been 
consistent with the outcome of the zero temperature calculation. However, one might equally 
well start from the expression (p3"D, and take the limit T — > after the integration over u has 
been performed. One may wonder whether the two limits are the same, that is whether the 
operations of taking T — > and integrating over u commute. It is easy to see that they do not, 
unless we keep a non-zero fermion mass. Let us see this considering the expression (p3|) for 
n = and k\ = 0: 



njf(0) = -2a J du J ^dve~Tm m (28) 
x Tcoth^y e ~^^ 2 \2-^—w 2 - 1 | 

We observe that, if m is zero, the proper time u in the integrand appears only in the combination 
uT 2 /\eB\ (with the exception of coth it). This suggests that we perform the natural change of 
variable u — > u\eB\/T 2 . The above expression becomes: 



n 33 



0,m- 



=o(0) 



X 



Aa\eB\ 

^e-^ta+i) 3 ( 2 mr 2 (2Z + l) 2 - 1 



00 , , / de.B| 

du coth ^ 9 
V 1 z 



(29) 



and shows that IlQ 3 m=0 (0) would be proportional to 1/T in the limit T — > where we have 
coth(w|e-B|/T 2 ) ~ 1. Thus we can interchange the limit T — > and the integration over the 
proper time u to find the correct zero temperature limit only if m 7^ 0, at least as long as 
\eB\ 7^ 0. This condition is consistent since the magnetic field always generates a dynamical 
mass when T < T c . 

Let us take the zero magnetic field limit of (^), after making the change of variable 
u — > \eB\u. We obtain then: 



/*00 P 1 2 

lim II 33 (A; ± ) = -2a / / c fo e - u [ ?n2+i ^( fc i+«'")] 

f (l-. 2 ) + i(2uWf + 



uvw n Wi 



(30) 



The integration in ( PU|) over the proper time u gives the relevant result of (equation (A12), 
after performing the momentum integration and changing the Feynman parameter x into (1 — 
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v)/2). For the case at hand {\eB\ = 0), we can take a massless fermion (m = 0) and the value 
of the photon thermal mass is found by setting first n = and then take the limit k± — > in 



M, 2 eB|=0 (r) = - Hmn^fcj.) = c aT 



k i ->0 



(31) 



with 



c = 4 



roc 

Jo ^ 



-u(2i+ly 



2(2Z + 1) 2 - 



u 



(32) 



To compute c, we use the Poisson resummation (^4j) and write 




= 8 In 2 (33) 

which gives the result that was found in || and [|14j] (with the notation a — > a /Air). 

Finally we decompose IT 33 (/c_i_) in a sum of two terms: the temperature independent part 
and the temperature dependent one. Using equation ( p3|) and the Poisson resummation (|24]), 
a straightforward computation leads to 

Uf(k ± ) = U° n (k ± ) + Ul(k ± ) (34) 
where IT°(fc_|_) is the zero temperature part ( p7|) and IT^(/c_|_) the temperature dependent part: 



—a 



7r eB J o 



f i 

/ r/r( T^bt .■, 



x 



X cosh it — cosh u^; u \ nv ,2 i 1 — Ij „ ,2l 
|efl| I" 1 i — W n\ 



k\ v " (cosh-uf — t> cothw sinhwf ) yY— lVe ^r 5 " cos[7mZ(l + u)]+ 



(35) 



sinh u 



i>i 



cothii 



\eB\ ^(-l) m e~i^ \ ~T2~ C0S [^ n K l + v )\ + 2nlitv sin[7rnZ(l + v)] 



We compute in appendix A the strong field (\eB\ — > oo) asymptotic form (|49|) of n° + n T that 
we used for the numerical analysis of (|). 

We note that the other components of the polarization tensor were computed |15| in 
3+1 dimensions and their computation in 2+1 dimensions would follow the same steps. 

Finally, we give the thermal photon mass, using the dimensionless variables already 
introduced: 



2 

^phot 



ton n ° 3 ^ 



2a r 00 



fc±^o \eB\ 



(36) 
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Figure 4: nl hot versus t for a = .01 and /i = .01 



where \i = m/y\eB\. We are actually interested in the behaviour of ITq 3 for all the values of 
the ratio \ijt: /x/t — ^ for the description of the phase transition where the dynamical mass 
vanishes for t = t c > and fi/t — > oo for the zero temperature limit and \i > 0. In figure § 
we plot fipfyjt as a function of t as well as its strong field asymptotic form obtained from the 
equation (|49|) of appendix A. 




t (1 + e-^) 2 ^ 

We can see the perfect agreement between the curves for the whole range of the ratio fi/t, as 
long as fi « 1 and t « 1 which is the case we study in the framework of the strong field 
approximation. We note the unusual behaviour of the thermal photon mass: for decreasing 
temperature y? phot increases as 1/t as long as t > /i, reaches a maximum and then decreses to 
when t « \i. The result is a large correction to the photon propagator when t ~ /i. 
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Figure 5: m^ yn ja versus T/a 



5 Phase diagram 

We wish to study the critical temperature and the critical field of the theory, defined by the 
relations m dyn (B,T c (B)) = and m dyn (B c (T),T) = 0. 

Let us first find the critical temperature when we fix the magnetic field. We plot in figure 
|5] the evolution of the dimensionless dynamical mass md vn /ot — Hn=o/& with the dimensionless 
temperature T/a. We start from T = and see that md yn first follows a plateau for small 
temperatures and then decreases when T reaches a value that we interpret as the critical 
temperature, T c . We cannot find any other solution than m dyn = when T > T c and the 
vertical slope when T = T c suggests that the transition might be of first order. We note that 
the recursion formula flT7| ) gives a second solution which is not physical since it does not lead to a 
converging iterative procedure in the integral equation ([)]) and that this unphysical solution also 
does not go beyond T = T c . Finally, we noticed that the convergence of the iterative procedure 
to solve the Schwinger-Dyson equation is slower as we approach the critical temperature. In 
the same figure we have plotted the results for two magnetic fields; the resulting dynamical 
mass and critical temperature are bigger for the biggest magnetic field, as expected. 

Let us now study the dynamical mass when we fix the temperature. We plot in figure 
H the evolution of m dyn /a as a function of yJ\eB\/a fixing the ratio T/a. We see that for a 
given non-zero temperature, there is a minimal value that the field must take to generate a 
dynamical mass, that is, a critical field. In the case T = 0, we see that no critical field is needed 
to generate a dynamical mass, as is known |§ (we didn't plot the curve down to \eB\ = since 
the strong magnetic field approximation is not valid in this region). Notice that the calculation 
of this curve, which has been done using the methods of this work (and taking the limit T — > 0) 
is identical to the corresponding result of the genuine T = calculation of ||. For T > 0, the 
transition again seems to be of the first order and the convergence of the iterative procedure is 
slower as we approach the critical field. If the temperature increases, the curve shows smaller 
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Figure 6: m^ yn ja versus \eB\^l 2 j 



a 



dynamical mass and higher critical magnetic field. 

We depict in figure [7| the phase diagram T c /a versus \J\eB c \/ot and make the com- 
parison with the constant mass approximation (where \i\ does not depend on I) discussed in 
appendix B. We see that the constant mass approximation systematically underestimates the 
critical temperature T C (B) and that this one is almost independent from the magnetic field for 
\eB\^l 2 jot > 100. Going beyond the constant mass approximation we find that the critical tem- 
perature increases with \eB\ and reaches asymptotically the value T c °° = a/2 when \eB\ — > oo. 
The fact that the temperature increases with increasing magnetic field is consistent with similar 
studies made in the Gross-Neveu model 1171 . 



It is interesting to compare the momentum-dependent case with the constant mass 
approximation along paths of constant temperature or constant magnetic field. In figure |8] 
we compare the constant mass approximation and the momentum dependent dynamical mass 
for a given value of the magnetic field. We see in the constant mass approximation that the 
dynamical mass is overestimated, as was found in ||, but that the critical temperature is 
smaller than in the momentum dependent study. 

In figure ^ we fix the temperature and plot the dynamical mass md yn /oi as a function of 
the magnetic field for both the momentum-dependent case and the constant mass approxima- 
tion. The constant mass approximation is seen to overestimate the dynamical mass for most 
of the range of the magnetic field; in addition the critical magnetic field turns out to be bigger 
for this approximation than for the momentum-dependent case. An important remark should 
be made here. For sufficiently high magnetic field strength the critical temperature is almost 
independent from the magnetic field. This means that if the temperature is set to a sufficiently 
large value, it will not be possible to find a corresponding critical field. Incidentally, this was 
the reason why we have chosen a safely low value for the temperature in figure ^. The results 
of figures [8] and [], are consistent with the phase diagram, where it was found that the constant 
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Figure 7: T c /a versus \eB\^l 2 ja 
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Figure 8: m dyn /a versus T/a for {eB^^/a = 100 
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Figure 9: rridyn/ a versus |e-B| 1//2 /a for T/a = 0.3. 



mass approximation underestimates the critical temperature (for constant magnetic field) and 
overestimates the critical magnetic field (for constant temperature). 



Finally we show in figure |T0| the ratio 2md yn (T = 0)/T c versus y\eB\/a, where md yn {T = 
0) is the dynamical mass that is obtained at the limit T = 0. This ratio is an important param- 
eter for superconductivity theories. We find out that this ratio is of order 10 and approximately 
constant over the whole range of y\eB\/ct considered here; this is consistent with the result 
obtained before in [[| where \eB\ = 0, which shows that the magnetic field has no important 
influence on this ratio. 



As a side remark, we computed the condensate 



0\W\o) = UmtrG(z,y) =txT^J T^f^M 



i - ( 27 

which is another order parameter for this transition. The LLL approximation gives: 



(38) 



4, /^ie-^ B lxffV 
J (2tt) 2 V 



M, 



Zfwf + Mf 



\eB\ 



71 



TJ2 



My 



(39) 



and we noticed numerically that < O^^O > /\eB\ almost does not depend on the temperature 
or the magnetic field, taking into account the momentum dependence of the self energy. When 
T > T c or \eB\ < \eB c \, the condensate is zero since the only solution for the series M; is 
Mi = for all I. The constant mass approximation is obtained very easily with the method of 
summation described in appendix B and the result reads 
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Figure 10: Ratio 2m^ n (T = 0)/T c versus \eB\ 1 / 2 / 



a 



(o|^|o) = ^tanh(^ (40) 

which shows that the condensate vanishes at the critical temperature or the critical field, 
when hq — 0. We note that an analytical computation of the condensate beyond the LLL 
approximation and the constant mass approximation has been performed in leading to 
(ffO|) when \eB\ oo. 



6 Discussion 

The effect of external magnetic fields on the state of the condensate is of great interest, in 
view of recent experiments with high-T c cuprates pertaining to the thermal conductivity of 
quasiparticle excitations in the superconducting state |0|. It is argued that quasiparticle 
thermal conductivity plateaux in a high-magnetic-field phase of these materials indicate the 
opening of a gap for strong magnetic fields, depending on the intensity of the external field, 
when the latter is applied perpendicularly to the cuprate planes. The phenomenon appears 
to be generic for systems of charged Dirac fermions in external magnetic fields, in the sense 
that strong enough magnetic fields are capable of inducing spontaneous formation of neutral 
condensates, whose magnitude scales with the magnetic field strength || 0, |20| . Such gaps 
disappear at critical temperatures of the order of the gap. 



The authors of |T8[ find that, for strong enough magnetic fields of 0[1] — O[10] Tesla, 
there are plateaux in the thermal conductivity of quasiparticle excitations about the (i-wave 
state (in particular about the nodes of the d-wave superconducting gap). Such plateaux are 
interpeted as an indication of the opening of a new gap, induced by the magnetic fields at the 
nodes. These plateaux disappear at a critical temperature that depends on the magnetic field 
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intensity, and in particular they report the empirical relation: 



T* xp oc yJ\eB\ (41) 

for the dependence of the observed critical temperature on the external magnetic field strength. 
The mass gap (that is, the dynamical mass at zero temperature) has been found in Q to 



be linear in y\eB\, so the critical temperature is expected to be linear in y \eB\. The above has 
been a conjecture from the zero temperature case, but a genuine finite temperature calculation 
had to be done before one could draw definite conclusions. The relevant results are contained 
in figure |7] and we may conclude the following: 

• In the constant mass approximation the critical temperature T c j a is constant (its value is 
about 0.4) for a wide range of the magnetic field. The critical temperature depends on the 
magnetic field only for relatively small field strengths. We note that this calculation has 
been done assuming a strong magnetic field, so one cannot explore the region of very small 
fields in a reliable way. Note that because of the choice of axes in figure |7| a square root 
dependence would show up as a straight line, which may be identified with the tangent 
to the curve. Thus the possibility that the behaviour of the critical temperature on the 
magnetic field conforms to the experimental result ( fill) is restricted to a small range of 



magnetic field strengths (approximately up to y\eB\/a ~ 50). 

In the case of the momentum- dependent fermionic self energy the situation changes quan- 
titatively: a straight line, corresponding to a square-root-like dependence fits to the data 



up to the value y\eB\/a ~ 150. Then the curve levels off and T c approaches 0.5a. which is 
somewhat bigger than the strong field limit in the constant mass approximation. However 
it is not clear what this "saturation" phenomenon implies for the superconductor case: if 
the field strength becomes very big, it will move the system out of the superconducting 
phase, so the whole calculation is not relevant any more. It appears that one should stay 
at a safe distance from both the "too weak" field region (where the LLL approximation 
collapses) and the "too strong" field regime, which the physical phenomenon may dis- 
appear Q For not too extreme fields the square root empirical relationship appears to 
persist in the momentum dependent case for relatively stronger fields than it is possible 
for the constant mass approximation. 

Acknowledgements 

This work has been done within the TMR project "Finite temperature phase transitions in 
particle Physics", EU contract number: FMRX-CT97-0122. The authors would like to ac- 
knowledge financial support from the TMR project. 

A Strong field approximation for the longitudinal polar- 
ization tensor 



Starting from the relations ( p7f ) and (35), we compute here the asymptotic forms of II and 
II T in the strong field limit where the dominant contribution in the integration over u comes 

4 We thank the referee for pointing this out. 



17 



from the region u » 1, which can be seen with the change of variable u — > u\eB\. We will use 
the dimensionless parameters previously introduced. For the temperature independent part, 
we can take fi = and we write that for u » 1: 



coshwf - v cothwsinhwt; - u (x-v) , /, \ -ud+v) 

~ (l-v)e {X ) + (l + v)e u ^ 1+v > 

sinhw 

k 2 k 2 

_ J_ cosh u — cosh uv _ _\_ 

g JeB\ 2sinhu rvJ g 2|eB (42) 



such that we obtain from 



-a: 



OO f\ 



k\e 2 i eS i / du dvy/ue u 4 



k\e 2^sT / dv- 



1 I 1 

1 



3/2 



1 - v)e-< x ~ v) + (l + v)e- 



-u(l+v) 



(43) 



For the temperature dependent part, we have to keep \i 7^ for the reason previously explained 
and we will consider the dominant term in the strong field approximation, i.e. 



Ti^lk I 1 — O < , 7,-1 r, r°° /"l fc ± cosh it -cosh ui> / 2 1 1-v 2 2\ 

Ylt - "AD" 1 ) 1 du dve-vm^^; u( M+ —u, n) (44) 

\&-t>\ \T^t JO ■/ — 1 



cothw 



-e 4 



! 2 

4^ cos[7rn/(l + d)] 



M 3/2 

eSl^-iy+H 2 f 1 dvcosfrrdtl+v)] f°° ^-e'^^^'^ 

J-l Jo U ' 



We first have the following integration over u which leads to the Bessel function K1/2: 

= f = 2 © V4 ^(2^) = f b e~^ (45) 

with a = fi 2 + ^-f-u) 2 and b — -rL. Then we approximate the integration over v by interpolating 
between the asymptotic expression of the integral when n » 1 and the one with n = 0. When 
n » 1, the oscillations of the integrand make the contributions of opposite sign cancel, such 
that the dominant contribution for the integral comes from the region 1 — < v < 1 where 
the cosine does not vanish: 



-1 



dv cos[7rn/(l + v)]e 



LJ.,2 , i=«f „,2 



cfc cos[7rnZ(l + f )]e~ T 

- _ 2HT 



2 ^ 

e * 



7r|n|/ 



(46) 
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We also have Jo(//, i) = 2e t and therefore we make the ansatz: 



2e 



I u 



1 + 7r|n|/ 

With these approximations, we can finally write 

\eB\ ~ t §f j l+7r\n\l 

The strong field asymptotic longitudinal polarization tensor is then, from (|4*3| ) and (f|S] 



(47) 



(46 



\eB\ 



-2ae 2 i eS i 



y/2- 



k 2 



eB\2 + u 2 



1 +fE(-i) m T ; 

1 1>1 i 



+ %\n\l 



e~ * 



(49) 



B Constant mass approximation 

We compute here explicitly the summation over the Matsubara modes in the integral equation 
(0) for the dynamical mass in the approximation where the latter does not depend on the 
Matsubara index I : \i\ = /x . We suppose that the correction to the photon propagator also 
does not depend on the Matsubara index and keep its transverse momentum dependence only. 
The integral equation (|9f) for the dynamical mass with the help of (|T3|) yields (we take Z\ — 1): 



/•oo „ 

1 = 2a / rdre~ r /2 
Jo 



Z](r, no) 



r 2 + fi 2 (r 



(50) 



where fj*(r) = r^hllf (r 2 \eB\) and the sums over Matsubara modes are 



St(r,/x ) = tJ2 



1 



S 2 (r,/i ) = *X! 



f (£ 2 + /^[(^-£ ) 2 +r 2 ] 
1 



7- (£? + //§)[(£, -£ ) 2 + r 2 F 



(51) 



has been computed in p2j and we shortly repeat here the steps of the computation. Ej can 
be expressed in terms of the contour integral 



If 1 

s t0^o) = ^— / dw— ; — 
2m Jc 1 + e 



1 



The contour C runs around the poles i(2l + \)xt of the Fermi-Dirac distribution function with 
residues equal to 1/t. This contour can be deformed to run around the 4 poles ±/xo and ±r+iu) 
of the rational fraction. One must not forget that this deformation leads to an overall minus 
sign since we travel along the new contour in the clockwise direction. The summation over the 
4 corresponding residues gives then 



1 . h (lM>P 
tanh — — 



(bl + r 2 -fi 2 



2/i C 



{u 2 + r 2 - fi 2 ) 2 + Afi 2 u 2 
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2r I 2 J (lu 2 -r 2 + fif + 4r% 2 



(53) 



where /3 = p^\eB\ = 1/t. 

£ 2 is simply given by S 2 (r, /i ) 



■ii^fo/V), so: 



Mo) 



1 i / /W 
tann 

2/io V 2 , 



[(r 2 - /i 2 + cu 2 ) 2 + 4/4^ 



2,:,212 



- r2 + Mo 



2 

+ ^coth *Z-r' + l4 

^4r 3 ^2 J (u 2 - r 2 + /i 2 ) 2 + 4r 2 ^ 2 

, 1 h (rj\ + /xg) ~ fog ~ r 2 + ^g) 2 

^ 2r ^ 2 y [{ul - r 2 + /i 2 ) 2 + 4r 2 ^ 2 ] 2 

The critical temperature T c in the constant mass approximation satisfies the equation 

r 2 /4(r) 



8r 2 sinh 2 (f ) foo - r 2 + /ig) 2 + 4r 2 ^ 2 

2 



with 



l=2a rdre- r /2 
Jo 



SL(r,0) 



r 2 + l^l(r) 



S 2 c (r,0) 



(54) 



(55) 



Si(r,0) 
S?(r,0) 



1 1 , /VflA cDg - r 2 

4 r 2 + c^ 2 2r C ° I 2 i {u 2 + r 2 ) 2 



/3 C 



4 (r 2 + cD, 



2A2 



+ 



Pc 



Cj 2 c — r 2 



8r 2 sinh 2 (^) fo 2 + 



,2\2 



J_ (rpA u 2 - r 2 J_ 3cDg - r 2 

+ 4r 3 C ° th I 2 J fo c 2 + r 2 ) 2 + 2r C ° th I 2 J (cu 2 + r 2 ) 3 



where Cj c = 7it c and fi 2 c (r) is the opposite of the strong field asymptotic form of the dimensionless 
longitudinal polarization tensor (|49| ) for zero Matsubara index taken at the critical temperature: 



Mr) 



ae 



ae 



-r 2 /2 



-r 2 /2 



^2r 2 + -limYi-lY+He- 1 ^ 
t c m-o ~ 



v/2r 2 + - 



(56) 



where the summation over I has been done to obtain equation (0). The solution of fl5"5|) is 
plotted in figure 0. 
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